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Local weather factors are widely considered to influence the transmission of infectious gastroenteritis. Few 
studies, however, have examined the non-stationary relationships between global climatic factors and 
transmission of infectious gastroenteritis. We analyzed monthly data for cases of infectious gastroenteritis 
in Fukuoka, Japan from 2000 to 2012 using cross-wavelet coherency analysis to assess the pattern of 
associations between indices for the Indian Ocean Dipole (lOD) and El Nino Southern Oscillation (ENSO). 
Infectious gastroenteritis cases were non-stationary and significantly associated with the lOD and ENSO 
(Multivariate ENSO Index [MEI], Nino 1 + 2, Niiio 3, Nino 4, and Nino 3.4) for a period of approximately 1 
to 2 years. This association was non-stationary and appeared to have a major influence on the synchrony of 
infectious gastroenteritis transmission. Our results suggest that non-stationary patterns of association 
between global climate factors and incidence of infectious gastroenteritis should be considered when 
developing early warning systems for epidemics of infectious gastroenteritis. 



Infectious gastroenteritis contributes significantly to the 1 billion episodes of diarrhea and 3 million deaths in 
children under 5 years of age per year, and is the fifth-leading cause of death worldwide' The transmission of 
infectious gastroenteritis is complex and multifactorial, involving both host and environmental factors. 
Local weather factors such as temperature, relative humidity, and rainfall have been suggested as important 
factors in the spread and seasonality of infectious gastroenteritis^ In addition to local weather factors, several 
studies reported that the El Niiio Southern Oscillation (ENSO) and Indian Ocean Dipole (lOD) play important 
roles in the transmission of infectious diseases, including dengue""", malaria'^"^' and cholera^^"^^. The ENSO is 
the most prominent source of interannual global climate variability which affects weather conditions, such as 
temperature, rainfall, wind speed and direction, and storm tracking throughout the world^"'^''. These effects, 
however, fluctuate and vary from region to region^"'^''. The lOD is another global climate phenomenon that arises 
from ocean-atmosphere interactions which affect climate patterns in the tropical Indian Ocean"'^"'^"''^*. Moreover, 
the World Health Organization (WHO) quantified the impact of global warming on diarrhea, and reported that 
warming by 1°C was associated with a 5% increase in diarrhea^^. Although regional differences and contrasting 
effects of temperature on different kinds of diarrhea are evident^", few studies have examined the non-stationary 
associations of infectious gastroenteritis and global climate variability. 

Wavelet analysis is useful in the investigation of non-stationary associations using time series data^' as it can 
measure associations (coherency) between two time-series at any frequency (period) band and time-window 
period. This analysis has been used to determine whether the presence of a particular periodic cycle at a given time 
in disease incidence corresponds to the presence of the same periodical cycle at the same time in an exposure 
covariate^". Wavelet analyses have been used to analyze the transmission of infectious diseases'''"''^'^^'^"'^'. 
Therefore a better understanding of the sensitivity of these analyses to climate variability might help develop a 
reliable climate-based prediction system for epidemics of gastroenteritis. 

Here, we explored the time-varying relationship between climate variation and monthly incidence of infectious 
gastroenteritis between 2000 and 2012 in Fukuoka, Japan. To our knowledge, this is the first report to quantify the 
time-varying impact of climatic factors on the number of infectious gastroenteritis cases using cross-wavelet 
analysis. 
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Results 

A total of 654,254 cases of infectious gastroenteritis from 2000 to 
2012 were included in our analyses, of which 392,514 (60.0%) were 
children aged 0 to 4 years, 171,750 (26.3%) aged 5 to 9 years, 47,541 
(7.3%) aged 10 to 14 years, and 42,449 (6.5%) aged 15 years or older. 

The time series for the number of infectious gastroenteritis cases 
per month, ambient temperature, relative humidity and rainfall dur- 
ing the study period are shown in Figure 1. As noted above, the 
incidence of infectious gastroenteritis displays a seasonal pattern in 
temperate areas, with marked peaks in winter (Fig. 1). 

The time series for DMI and ENSO indices (MEI, Nifio 1 + 2, 
Nino 3, Nino 4, and Niiio 3.4) during the same period are shown in 
Figure 2. Strong positive lOD events (indicated by large DMI values) 
occurred in 2006, with a peak DMI in October, and in 2012, with a 
peak DMI in August. Strong ENSO events (indicated by large MEI 
values) were observed in 2006 and again in 2009 to 2010 (Fig. 2). 

Cross-wavelet coherence and cross-wavelet phase analysis of the 
global climatic time series (DMI and ENSO indices) with infectious 
gastroenteritis cases by month are shown in Figure 3. The wavelet 
coherence provides information about whether two non- stationary 



signals are linearly correlated at a particular time and frequency"'^". 
The wavelet coherence is equal to 1 when there is a perfect linear 
relationship at a particular time and frequency between the two 
signals'"*'". In fact, cross-wavelet coherency analysis revealed that 
infectious gastroenteritis cases were significantly (p < 0.05) coherent 
with DMI for 2 years (2006-2007) and 1 year (2010). With regard to 
ENSO indices, MEI was significantly coherent at periods of approxi- 
mately 1 to 2 years (2005-2006). Niiio 1 + 2 was significantly coher- 
ent at a period of approximately 1 year (2003-2004) and 
approximately 2 years (2006). Niiio 3, Niiio 4, and Niiio 3.4 were 
significantly coherent at a period of approximately 1 to 2 years 
between 2003 to 2008. DMI and ENSO indices were predictive of 
disease data, as shown by the cross-wavelet phases (Fig. 3). AU local 
weather factors had non-stationary relationships with both DMI and 
ENSO indices (Supplementary Fig. SI). 

Discussion 

In this study, we observed that the monthly incidence of infectious 
gastroenteritis cases is significantly associated with changes in the 
lOD and ENSO (MEI, Nifio 1 + 2, Niiio 3, Nifio 4, and Niiio 3.4) for 
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Figure 1 | Monthly time series data in Fukuoka, Japan (2000-2012). (a) Infectious gastroenteritis cases; (b) Temperature; (c) Relative humidity; (d) 
RainfaU. 
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Figure 2 | Monthly time series data for global climatic indices (2000-2012). (a) DMI (dipole mode index); (b) MEI (multivariate ENSO index); (c) Nino 
1 + 2 (ENSO index); (d) Niiio 3 (ENSO index); (e) Nino 4 (ENSO index); (f) Niiio 3.4 (ENSO index). 



the 1- to 2-year periodic mode. This association is non- stationary 
and appears to have a major influence on the synchrony of infectious 
gastroenteritis transmission. These findings indicate that early warn- 
ing systems for epidemics of infectious gastroenteritis should con- 
sider non-stationary, and possibly non-linear patterns of association 
between climatic factors and infectious gastroenteritis cases. 

Predicting the potential effects of changes in climate on the incid- 
ence and distribution of infectious gastroenteritis requires an under- 
standing of the relationship between non-stationary climate and 
transmission of disease. Our results demonstrate that the incidence 
of infectious gastroenteritis is strongly associated with non- station- 
ary climate variables, including DMI and ENSO indices, with coher- 
ent cycles of approximately 1 to 2 years. A study in Peru reported that 
during El Niiio there was a marked increase in the number of cases of 
diarrhea and dehydration in infants and young children'"'. Further, 
another study in Peru indicated that the number of hospital admis- 
sions for infectious gastroenteritis increased during the El Niiio 
event, and that this increase was greater than that explained by 



seasonal temperature changes^. Thus, our results are consistent with 
previous studies and demonstrate the importance of climate change 
in the spread of infectious gastroenteritis. 

Extreme lOD and El Niiio are two dominant drivers of year-to- 
year climate variability''. In recent observations, lOD and El Nino 
repeatedly co-occurred, indicating the interactive nature of the two 
major climate modes''*. A recent study indicated that the warming of 
the Indian Ocean relative to that of the Pacific might play an import- 
ant role in the climate changes of the latter". Another study sug- 
gested that the Indian Ocean might significantly enhance El Niiio 
and its forecast of onset, which in turn enhances the lOD and its 
long-range predictability". Thus, our combined lOD and ENSO 
results demonstrate the importance of global climate factors on the 
transmission of infectious gastroenteritis. 

Our study also found that the incidence of infectious gastro- 
enteritis is strongly associated with local weather factors, such as 
temperature, relative humidity, and rainfall, with coherent cycles 
throughout the year. The positive relationship between infectious 



SCIENTIFIC REPORTS | 4:5157 | DDI: 1 0. 1 03 8/srep051 57 



3 



(a) 



Coherence 




Time 
Phase 




Coherence 




I 



Time 
Phase 




I 



Time 



Time 



(c) 



Coherence 




Time 
Phase 




Coherence 




I 



Time 
Phase 




I 



Time 



Time 



(e) 



Coherence 




Time 
Phase 




(f) 



Coherence 



I 



Time 
Phase 




I 



Time 



Time 



Figure 3 | Cross-wavelet coherency and phase of the infectious gastroenteritis time series with global climatic indices, (a) DMI (dipole mode index); 
(b) MEI (multivariate ENSO index); (c) Nino 1 + 2 (ENSO index); (d) Nino 3 (ENSO index); (e) Niiio 4 (ENSO index); (f) Niiio 3.4 (ENSO index). Red 
regions in the upper part of the plots indicate frequencies and times for which the two series share variability. Cone of influence (within which results are 
not influenced by the edges of the data) and significant coherent time-frequency regions (p < 0.05) are indicated by solid lines. In cross- wavelet 
phase plots, colors correspond to different lags between the variability in the two series for a given time and frequency, measured in angles from —PI to PI. 
A value of PI corresponds to a lag of 16 months. 
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gastroenteritis cases and temperature, humidity, and rainfall in our 
present study is broadly consistent with previous studies in Japan, 
China, Peru, Fiji, and Bangladesh^*' " '^"'''. Further, temperature 
increases were suggested to be positively associated with Salmone- 
lla, Campylobacter, and Escherichia co/P*'''^ Relative humidity and 
rainfall also affect water and sanitation infrastructure and the num- 
ber of pathogens, and might impact the replication rate of certain 
bacterial and viral pathogens^". Our results suggest that previous 
findings based on this assumption can be improved by assuming a 
non-stationary association and more accurately evaluating the pos- 
sible non-linear association between climatic factors and infectious 
gastroenteritis transmission. 

Although data on the cause of infectious gastroenteritis were not 
available for analysis, a striking increase was suggested for both the 
number of cases of bacterial diarrhea during summer and of rota- 
virus during winter"". Hot weather might cause certain behavioral 
patterns, such as higher water consumption and less conscientious 
hygiene practices, which promote diarrhea transmission"". In addi- 
tion, this study found long-term lag effects of global climate variabil- 
ity on infectious gastroenteritis transmission, which might be related 
to the fact that hot weather affects contamination in food production 
and distribution systems"'""'^. This suggests that temperature control 
during food production, processing, transport, preparation or stor- 
age might interact with weather patterns and thus have an impact on 
the risk of disease. For childhood infectious diseases, another study 
also suggested that periodic diseases become chaotic with introduc- 
tion of low-level vaccination"*^. The Ministry of Health, Labour and 
Welfare, Japan has not adopted routine rotavirus vaccination, and 
the results of this study might be affected by the effects of rotavirus 
vaccination. This characteristic may also indicate the need for more 
precise modeling to elucidate the complex relationships of global 
climate variability and infectious gastroenteritis transmission, and 
further investigation would be a critical topic for future study. 
Although other behavioral factors unrelated to changes in climate 
might also affect the seasonality of diarrhea, such as changes in 
patterns of food availability, we have no evidence for such factors. 

Our findings with the lOD and ENSO indices support the recent 
proposal that large-scale climate indices might be more useful for 
forecasting ecological effects than local weather variables"*^. Climate 
can affect the transmission of infectious diseases through several 
linear and non-linear pathways in a host population. Climate can 
also affect several biological traits of the organisms involved in the life 
cycle of parasites, from individual life histories to population 
dynamics"*"*. The WHO and the Intergovernmental Panel on 
Climate Change (IPCC) have identified changes in the incidence of 
diarrhea as one of the most important future health effects of climate 
change"*^. While local weather most likely only affects the biological 
components of disease transmission, large-scale climate patterns 
might also influence contextual components of disease transmission, 
such as population susceptibility. 

Our results for infectious gastroenteritis in this paper do not con- 
sider etiology despite the fact that various pathogens have different 
seasonal patterns. As the response of bacterial and viral pathogens to 
the effects of climate variability might differ, our results might only 
reflect the responses of the most common cause of infectious gastro- 
enteritis. Our results should therefore be interpreted with caution, 
and further studies on the effects of climate variability on specific 
pathogen-induced infectious gastroenteritis are warranted. 

The present study has several limitations. First, not all cases in the 
community are represented in the surveillance data. This under- 
reporting can occur anywhere in the reporting chain, from the initial 
decision of a patient not to seek health care to failure to record cases 
in the disease registry. However, we consider this did not result in 
substantial bias because the degree of under- reporting is unlikely to 
have varied over time. Second, the participating sentinel medical 
institutions were recruited on a voluntary basis. However, this did 



not pose a threat to the validity of the comparisons over time, which 
is the subject of this study. Third, we analyzed monthly data, span- 
ning from 2000 to 2012. Modeling accuracy of these associations 
would be improved by a longer study period or more detailed data. 

The results of this study have practical implications for public 
health officials. Elucidation of the relationship between climate vari- 
ability and infectious gastroenteritis transmission is important for 
disease control and prevention. The present findings may help public 
health officials to predict epidemics and prepare for the effects of 
climatic change on infectious gastroenteritis through the imple- 
mentation of preventive public health measures, such as warning 
health workers, active community-based campaigns promoting oral 
rehydration, increasing staff and the number of beds in oral rehydra- 
tion units in hospitals, early warning weather forecasts, and improv- 
ing infectious gastroenteritis control programs. 

In conclusion, this study found quantitative evidence that non- 
stationary associations have a major influence on synchrony between 
changes in both the lOD and ENSO and the monthly incidence of 
infectious gastroenteritis cases for the 1- to 2-year periodic mode. 
Our results suggest that non- stationary patterns of association 
between global climate factors and infectious gastroenteritis cases 
should be considered when developing early warning systems for 
epidemics of infectious gastroenteritis. 

Methods 

Data sources. Under the Infectious Disease Control Law, the number of infectious 
gastroenteritis patients is reported on a weekly basis from 120 sentinel medical 
institutions in Fukuoka Prefecture in southwest Japan^^'^^. A case of infectious 
gastroenteritis is defined by clinical factors of sudden stomach ache, vomiting, and 
diarrhea. However, under this law it is not necessary for sentmel medical institutions 
to report information from laboratory tests. We therefore obtained clinical data that 
were recorded and reported by sentinel volunteers to the Fukuoka Institute of Health 
and Environmental Sciences, which is the municipal public health institute of the 
Fukuoka Prefectural Government. JMonthly cases of infectious gastroenteritis were 
calculated from daily records based on the day of diagnosis. Large outbreaks were 
removed from the data set. 

The strength of the lOD was measured by the Dipole JVIode Index {DJMI), defined 
as the difference in sea surface temperature (SST) anomalies between the western 
(lO'S-lO'N, 50 -70 E) and eastern (lO^S-O', 90°-110°E) tropical Indian Oceans"". 
DMI data were obtained from the Japan Agency for Marine-Earth Science and 
Technology (JAMSTEC) (http://www.jamstec.go.jp/frcgc/research/dl/iod/). The 
ENSO indices (i.e. JVIultivariate ENSO Index [MEI], Niiio 1 + 2, Niiio 3, Niiio 4 and 
Niiio 3.4) were also extracted from the National Oceanic and Atmospheric 
Administration (NOAA) Climate Prediction Center data (http://www.cpc.ncep.noaa. 
gov). 

We also obtained data on daily average temperature, relative humidity and rainfall 
in Fukuoka Prefecture from the Japan JMeteorological Agency. iVIonthly means for 
average temperatures, relative humidity and rainfall were calculated from daily 
records. 

This study was approved by the ethics committee of the Fukuoka Prefectural 
Government. The requirement for written informed consent was waived. Patient 
records and other patient information remained anonymous and de-identified prior 
to analysis. All methods were carried out in accordance with approved guidelines and 
regulations. 

Statistical analysis. We examined the frequency of reported infectious gastroenteritis 
cases using cross-wavelet coherency analysis^^'^^--*^-™. This analysis enables the 
investigation and quantification of whether the presence of infectious gastroenteritis 
at a particular frequency at a given time corresponds to that of a climate covariate at 
the same frequency and time. Specifically, cross wavelet coherency analysis was used 
to determine whether a particular frequency of reported cases at a particular time 
during an infectious gastroenteritis epidemic corresponded to the same frequency of a 
climate covariate at the same time. Coherence analysis is similar to correlation 
analysis as it is normalized between 0 and 1, where 0 corresponds to the total absence 
of cycles with the same period in the analyzed time series, and 1 corresponds to the 
presence of cycles with exactly the same periods in the analyzed time series^^'^^-^^-'"'-'*'. 
Briefly, for the two time series, the cross wavelet spectrum (CWS) is estimated over 
time for a series of frequencies (or period, i.e., 1/frequency). The CWS is normalized 
by the power spectra of two variables in the same manner that covariance is 
normalized by the product of the standard deviation of two variables when their 
correlation is estimated. Therefore, CWS can be considered equivalent to covariance. 
When the magnitude of the CWS is similar to the normalized product of the time 
series spectra the ratio is equal to 1, which indicates that the two studied time series 
have cycles with exactly the same periodicity*^'"*^'*^^'"*"'*^. 

We also estimated the cone of influence, where inferences from the wavelet analysis 
outside the cone are not valid because of the manipulations performed to generate the 
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wavelet spectrum in the absence of data for more frequent periods Cross 
wavelet coherence significance was estimated using the method described by Maraun 
and Kurths*^ for a minimum of 6 months (i.e., the minimum period of interest in 
cycles studied with cross wavelet analysis was 6 months), with a total smoothing 
window of 31 (i.e., 15 months) and dimensionless parameter of 6 for the Morlet 
Wavelet. Further details concerning cross wavelets and software are described by 
Maraun and Kurths*^ and Chaves and PascuaP". All statistical analyses were con- 
ducted using R version 3.0.2 (R Core Team, R Foundation for Statistical Computing, 
Vienna, Austria). 
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